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~j ' Abstract 

Ch ■ We construct a Poisson map between manifolds with linear Poisson brackets corre- 

sponding to the two samples of Lie algebra e(3). Using this map we establish equiv- 
alence of the Steklov-Lyapunov system and the motion of a particle on the surface of 
the sphere under the influence of the fourth order potential. To study separation of 
variables for the Steklov case on the Lie algebra so(4) we use the twisted Poisson map 
between the bi-Hamiltonian manifolds e(3) and so(4). 
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1 Introduction 

A standard form of the Kirchhoff equations is the following: 

M = Mxn + pxu, p = pxn, (1.1) 

where X x Y stands for the vector product of three-dimensional vectors jH]- 

These equations describe the motion of a rigid body in the ideal incompressible fluid if 

two vectors M and p are the impulsive momentum and the impulsive force while the vectors 

fl and U are the angular and linear velocities of the body. All these vectors in M^ are 

expressed in the body frame attached to the body [HI. 

It is known [llj that the system ()1.1|) is Hamiltonian with respect to the Lie-Poisson 

bracket of the Lie algebra e(3) of the Lie group E{3) of Euclidean motions of M^ , i.e. with 

respect to the Poisson bracket 

{Mi, Mj}i = EijkMk, {Mi,pj}i = EijkPk, {Pi,Pj}i = , (1.2) 

where Sijk is the totally skew-symmetric tensor. Here and below we identify dual space e*(3) 
with e(3) using the standard inner product [TT] . 

The Hamiltonian equations of motion for an arbitrary Hamilton function H = H{p, M) 
in the bracket ()1.2|) read 

M = M X VmH + px VpH, p = px VmH. (1.3) 

This Euler's equations on e*(3) coincides with ()1.1|) if 

H{p, M) = (M, AiM) + (M, A2J9) + (p, A3P). 

Here A^ are special numerical matrices and (., .) stands for the standard scalar product in 

R^ HI]. 

The Poisson bracket ()1.2|) has two Casimir function 

A=\p\\ B={p,M). (1.4) 

which are therefore integrals of motion for p.lj) in involution with Hamiltonian H{p, M) and 
with any other function on the phase space. 

The Steklov-Lyapunov case of the rigid body motion (the Steklov-Lyapunov system, for 
brevity), is characterized by the following diagonal matrices 

A = diag(ai,a2,a3), C = diag(a2 - 03, ^3 - ^i, cti - ct2) ■ (1-5) 

In pi] Steklov found a Hamilton function 

4Hs{p,M) = (M,AM) + 2(M,A^p) + (Ap,C2p) (1.6) 

for which equations (jl.3p possess a fourth additional integral. Here wedge denotes adjoint 
matrix, i.e. cofactor matrix. In our case it reads A^ = (det A) A~^. 

Later Lyapunov independently discovered an integrable case of the Kirchhoff equa- 
tions whose Hamiltonian 

4Hl{p, M) = (M, M) - 2(M, Ap) + {p,C^p) (1.7) 



is a linear combination of the Steklov integrals ^3] and the Casimir functions (jl.4p . 

The Lax matrices for the Steklov-Lyapunov system may be extracted from the Kotter 
work m 12] . Namely, in [7] Kotter introduce two vectors I and m 

£(A) = W(^(M-Bp) + Aj9j , m(A) = W^p, (1.8) 

where W and B are diagonal matrices with the following entries 

n=3 

Wji = a/a -a,, Bji = ^ l^ij-fcl ttfc . (1.9) 



Here A is auxiliary variable (spectral parameter) and functions Wjj = v^A — a^ can be con- 
sidered as basic elliptic functions (see [I1I21)- The equations of motion 

^HX) = {Hl, eh = miX) X ^(A) , (1.10) 

at 

may be rewritten in the Lax form 

^^e(A) = [^e(A),^e(A)] (1.11) 

using two Lax matrices 

3 3 

^e(A) = 5^i'.(A)a, , ^e(A) = J2m,{X)ai , (1.12) 

where cjj are the Pauli matrices. 

For an actual integration of the corresponding Hamiltonian flow in terms of elliptic 
functions see [7J, and for a more modern account P3I21- In fact in the known integration 
procedure we have to use the Lax matrices ()1.12|) with the spectral parameter A varying on 
an elliptic curve. 

However, it is known that the Steklov-Lyapunov flow is linearizable on the Jacobian of 
hyperelliptic curve instead of elliptic curve, similar to the Neumann system 0. Therefore, we 
can suppose that the Steklov-Lyapunov system belong to the family of the Stackel systems 
Pnj and there are the Lax matrices with rational dependence on the spectral parameter. 

The desired rational Lax matrices were constructed by Bolsinov and Fedorov by using 
a special triplets of vectors, which are coordinates on some comphcated phase space (see 
and references within). 

The aim of this note is to identify the Steklov-Lyapunov system with an integrable 
motion of a particle on the surface of the sphere, which is the Stackel system. This result is 
the direct sequence of the Kotter separation of variables j7] . 

2 Separation of variables 

If (Tfc are 3x3 Pauli matrices, then the Lax matrix =Sfe(A) (|1.12|) looks like 

4 -£2 
^e(A) =1-40 ii \ e so{3) . 
£2 -ii 



The corresponding spectral curve is defined by equation det(/i— =Sfe(A)) = 0, which is reduced 
to the following equation 

C: ^2^P3(A) = 0, (2.1) 

where 

F3(A) = |£|2 = a^\^ + (/3 - ahi X)X^ + H1X + H2. (2.2) 

Coefficients Hi^2 are linear combination of the Steklov- Lyapunov integrals ()1.6II1.7|1 and the 
Casimir functions ()1.4|1 

B B 

Hi = Hl trA + ^trA"^, i/2 = -^5 + - tr A^ - ^ det A . (2.3) 

The separation of variables associated with the curve ()2.H) was constructed by Kotter 
[7j and may be recovered in framework of the modern Sklyanin method |13j. 

Proposition 1 Separated variables Ui^2 associated with the Lax matrix =^e(A) are poles of 
the corresponding Baker-Akhiezer function ^(A) with the following dynamical normalization 

" = T7V-r — ^-^ r Wp, a G R . (2.4) 

|(M-Bj9)xp| ^' ^ ^ 

The proof consists of direct comparison of the known separated variables ^ with poles of 
the Baker-Akhiezer function ^(A), which is an eigenvector of the Lax matrix 

^e(A)*(A) = /i*(A). (2.5) 

Since an eigenvector is defined up to a scalar factor one has to fix a normalization of ^(A) 
imposing a linear constraint 

(a,*(A)) = l. (2.6) 

In our case =Sfe(A) G so(3) and excluding /x from ()2.5ll2.(ij) we derive that poles Uk of ^(A) 
are roots of the following equation 



.a. X l,Q. X 



')=0, 



where i is the Kotter vector (jl.8|) . 

Inserting a (12. 4j) in this equations and dividing it on polynomial det(A — A) one can 
define the separated variables ui^2 as roots of the following function 

^^^^ _ {X-ui){\-U2) _ {axe,axi) _^ x\ _^ 

^^~ det(A-A) " det(A-A) " Z^ A - a, " "^ ' ^ '^ 



where vector x is given by 



_ (M - Bp) X p 

\{M -Bp) X pI ' 



The equation ()2.7|) coincides with definition of the separated variables from [7] , where Kotter 
also proved that initial equations of motion (|1.1|) can be written in the form 

2v/P3(«i)det(A-Mi) . 2y/Ps{u2) det(A - U2) ,^ ^, 

Ml = , U2 = (2.8) 

Ui -U2 Ui - U2 

and then he integrated these equations using Abel-Jacobi inversion theorem. 



^.= -^{^i,e(A)} 



Proposition 2 The Kotter variables Ui2 ji^. 7| ) and momenta Vi^2 defined by 

J = 1,2, (2.9) 

are canonical variables 

{ui,U2}l = {^^1,^2)1 = 0, {Vj,Uk}l = Sjk, 

which satisfy to the following separated equations 

The proof is straightforward. 

As sequence of ()2.10p . in canonical variables Ui^2 and fi^2 integrals of motion Hi 2 are the 
Stiickel integrals 

2 

Hj = E '^7^' (^(^'^■) ^^ - ^(^^)) ' ^' = 1' 2' (2-11) 

A:=l 

where 

(/?(A) = det(A - A) , U{X) = a^X^ + {/3 - ahr A)X^ , 

and S is the Stackel matrix 

5= (7 "^). (2.12) 

This matrix is 2 x 2 block of the transpose Brill-Noether matrix Uc 

which determines the Abel- Jacobi map on Jacobian of hyperelliptic curve C fl2.ip completely. 

In our case the Stackel matrix S ()2.12|1 is the lowest block of the transpose Brill-Noether 
matrix Uc ()2.1Hj) ) and, therefore, there are canonical coordinates in which equations of 
motion ()2.8|) are the Newton equations [13 HZj . 

In the next section we discuss this canonical change of variables {p, M) -^ {x, J) which 
transforms the Kirchhoff equations p.lj) to the Newton equations in detail. 

3 The Poisson map 

The Poisson manifold is a smooth manifold Ai endowed with the Poisson brackets {. , .}m- 
If A4 and A/" are Poisson manifolds, a smooth map f : Ai ^ N" is called a Poisson map 
provided that it preserves Poisson brackets, i.e. 



for all 95, ■0 G C°°{J\f). Here f*ip = y? o / is a lifting of the function ip G C°°{Af) on M. 



Below we deal with linear Poisson brackets corresponding to the two samples of Lie 
algebra e(3) and to one sample of so(4) algebra, i.e. with homomorphisms of these Lie 
algebras. For brevity we will use the same notations both for the Poisson manifolds and the 
Lie algebras. 

If A4 coincides with A/", the Poisson maps are called canonical transformations. The 
problem of complete efficient description of all nonlinear canonical transformations is un- 
solvable. The reason is that for any function /(xi, . . . ,a;„) the flow defined by ODEs 
ii = {/ ) 2;i, . . . , Xn} yields a one-parameter group of canonical transformations. However 
one can investigate some interesting subgroups of nonlinear canonical transformations. 

Let X and J are coordinates on the Lie algebra e(3) with the standard Lie-Poisson 
brackets 

The brackets (j3.1|) respect two Casimir elements 

n=3 n=3 

A=\x\^ = Y,4, B={x,J) = Y,XkJk. (3.2) 

fc=l k=l 

Fixing values of the Casimir elements one gets a generic symplectic leaf of e(3) 

Oab = {x,J : A = a\ B = b}, (3.3) 

which is topologically equivalent to cotangent bundle T*S'^ of the sphere 

S^ = {x gMMx| = a}. 

Symplectic structure of Oab is different from the standard symplectic structure on T*S'^ by 
the magnetic term proportional to b JT] . 

If 6 = there is standard Poisson map T*S'^ -^ e(3) 

p : (tt, x) — i> J = vr X x, (3.4) 

where vr G M^ is conjugated to x momenta 

{7Vi,Xj} = Sij, and (a;,7r) = 0. (3.5) 

Let us consider classical counterpart of the Fourier transformation f : T*S'^ ^ T*S'^ 
defined by 

/: (x, tt) ^ (-TT, x) . (3.6) 

This symplectic mapping may be lifted to the Poisson mapping / : e(3) — > e(3) 



M xp M 

|M|2 |Mxp| ^ Jm] 



/: {p,M)^{x = ^la^-j^r^,,^ + bj—^, J = m]. (3.7) 



Its inverse mapping looks like 



/-'■■ (-'^)-(^=-^/«'-^^^ + /^^' M=J]. (3.8) 



The maps / and / ^ couple coordinates (x, J) on algebra e(3) (j3.ip with the following values 
of the Casiniir functions ()3.2j) 

A = a^ B = b 

and coordinates {p, M) on another sample of e(3) ()1.2|) with the following values of the 
Casimir elements (jl.41) 

A = a\ 13 = p. 

The symplectic map / fl3.6|) on T*S'^ is easy generalized: 

fg- ix,rc) ^ {-11 + g{x),x) 

if g{x) is a function on x such that 

{x,g{x)) = 0. 

For instance we can put g{x) = x x Bx, where B is an arbitrary numerical matrix. 

For brevity the lifting of this symplectic map fg to the Poisson map we present at the 
special case B = diag(6i, 62, ^a) and b = only. It will be enough in order to identify the 
Steklov-Lyapunov system with another integrable systems on e(3) |19j . 

Proposition 3 Let B and C are numerical diagonal matrices with the following entries 

n=3 

Bii = bi, Cu= ^ Sijk bj , keR. (3.9) 

j,k=i 



The mapping fg : {p, M) —>■ {x, J) defined by 

"'''" x,x X p 
is a Poisson map fg : e(3) — > e(3) such that 



(M-Bp)xp ^ ,^ 2^ 

x = a-i— — -^ — ^, J = M + a~^C 

M-Bp xp ' 



(3.10) 



A=\x\^ = a^, B = {x,J) =0. 
Here [y,z]^ is an " anticommutator" of two vectors y and z defined by 

n=3 
j,k=l 

The proof is straightforward verification of the Poisson brackets. 

Using B and C (j3.9p we determine matrix A = diag(ai, 02, as) and its adjoint A^ 

A = ^tr(B)-B and A"" = hc^ -B^) . (3.11) 

In these notations let us consider motion of a particle on the surface of the unit sphere 

52 = {xeM^ |x| = 1} (3.12) 
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under an influence of tlie fourtli order potential. The Hamilton function is equal to 

Hi{x,J) = j{J,J) -{x,Bx)(-f + d{x,Bx)) +S{x,A''x) , (3.13) 

where 7,^ are arbitrary parameters. If 5 = the function Hi{x,J) ()3.13p is Hamiltonian 
for the Neumann system. Recall that we can put |x| = 1 without loss of generality using 
canonical transformations x -^ a^^x. 

Proposition 4 Let o? and (3 are the values of the Casimir elements A and B \1.4^ on e(3) 
algebra M.^) . If 

7 = /3 - a^tr (A) , 6 = a^ (3.14) 

the Poisson map fg hS.lO^) identifies the Hamilton function Hi \S.1S\] on T*S'^ with the 
following Hamiltonian on the algebra e(3) ( li.iJ)) 



Hi = fg[Hi^ =i7z.-^/9trA + a2trA^ 



This Hamiltonian coincides with the Lyapunov integral H^ (|i.7| ) up to the last two terms 
depending on the Casimir elements. 

The proof is straightforward. 

As a result we arrive at the following conclusion: the Steklov-Lyapunov system is equiv- 
alent to the potential motion of a point a; G M^ constrained to the sphere, which belongs to 
the Stackel family of integrable systems jTHj. 

For the uniform Stackel systems we know the Lax matrices, the classical r-matrices, 
the bi-Hamiltonian description, the Backlund transformations, the separated variables, the 
theta-f unction solutions and many other facts. The Proposition 0] allows to transfer all these 
know results concerning to the Stackel systems onto the Steklov-Lyapunov system directly . 

4 The rational Lax matrices 

It is known [201, that the Hamiltonian Hi{x, J) ()3.13p is separable in elliptic coordinates Ui^2 
on the unit sphere S"^ ()3.12|) . which are roots of the standard form 

According to [1^ El ^] the generic 2x2 Lax matrices for the uniform Stackel system 
are constructed using Hamiltonian Hi and the generating function e(A) of the separated 
variables only 



^r{X) 



-}-eu{\) + w{X) e{\) ^e,(A) 




.(A) = . (4.16) 



Here et = {Hi,e} and function w{X) is given by w{X) = </)(A)e(A) ^ , where (f){X) is a 



MN 



parametric function and [^Ja/a^ is the hnear combinations of the following Laurent projections 

[HHII 



[e] 



Af 



E ^^^' 



A:=— oo 



Af 



E &^''- 



(4.17) 



J Af 



k=-M 



In our case N=2, M = and 0(A) = — 5A — 7. It's easy to prove that the required Lax 
matrices ()4.16p are given by 



n=3 / 1 ^^^^ 

if.(A) = J2 



i=l 



xf 



\ 



2 A — Oj A — ttj 

1 TTj^ 1 XjTTj 

V 4 A - a,- 2 A - a,- / 





(5A + (5(a;, Bx) + 7 



(4.18) 



where XjjTTj are canonical coordinates (j3.5p and 







S^rW 



(A - (x, Ax)) {5X + 5{x, Bx) + 7) - 6{{x, Axf - (x, A^x)) 



The Poisson bracket relations between the entries of the matrix J^f^iX) are closed into the 
standard linear r-matrix algebra. The corresponding r-matrix is rational dynamical matrix 
fnillZj- Here we present only one relation 



-- Ci A , e // \ = 

2 J^ X-n 



(4.19) 



which allows us to introduce canonical momenta 



Vk = - 2^*(^) 



X=Uk 



-, n=3 



X i 71 { 



2 ^ A - ai 

1=1 



X=Uk 



such that {t'fc,tfj} = (Jfcj- 

Inserting X = Uk into the determinant of the Lax matrix ^^(A) ()4.18|) 



det^,(A) 
one gets two equations 



1 2 1 
■-et + -eett-we' 



PsW 



det(A - A) 



_1 2.~ ^ _ Psjuk) 



k = l,2. 



where P3 is generating function of integrals of motion 

P3(A) = 5A=^ + 7A' + Ai?i + ^2. 
In Hamiltonian variables Vk, Uk these equations are the separated equations 



^ 



det{uk — A) 
9 



A: = l,2, 



(4.20) 



(4.21) 



whereas in Lagrangian variables u^, Uk they are equations 



1,2, 



det{uk — A) / det(«fc — A) 

which may be integrated using Abel-Jacobi inversion theorem. 
The second integral of motion H2 in ()4.2()|1 is equal to 

^2(x, J) = -i(J, AJ) + (x, A^x) (7 + S{x, Bx)) . (4.22) 

As functions on elliptic coordinates Ui^2 and momenta ■(;i^2 integrals of motion Hi 2 (I4.2U|) on 
T*S^ coincide with the Stackel integrals (j2.11|) . 

So, if conditions (|3.14|) hold then the Poisson map fg (J3.1U|) identifies the second integral 
of motion H2 fl4.22|) on T*S'^ with the following integral on the algebra e(3) ()1.2|) 

H2{p, M) = fg (#2) = -Hs + ^f3tT A^ - a' det A . 

This is the Steklov integral Hs fll.6|) up to the constant terms depending on the Casimir 
elements. 

The same Poisson map fg ()3.10|) identifies standard elliptic coordinates -Ui 2 (I4.15P on the 
sphere S"^ and the corresponding separated equations ()4.21|) with the separated variables Ui^2 
(j2.7p and the separated equations (j2.1Up for the Steklov-Lyapunov system proposed by Kotter 
[7]. It allows us to use the standard finite-band integration technique for the potential motion 
on the sphere ^^ in order to verify the Kotter solution of the Steklov-Lyapunov system in 
terms of the theta-functions. 

Applying the Poisson map fg ()3.10|1 to the Lax matrix ^^(A) ()4.18j) one gets the rational 
Lax matrix for the Steklov-Lyapunov system. In [Hj the similar Lax matrix was constructed 
by using the triplet of the vectors x, y, v defined by 

(M — 'Bp) = xxy, p = X X V . 

5 The Steklov system on so(4) and twisted Poisson map 

The bi-Hamiltonian manifold is a smooth manifold A4 endowed with a pair of compatible 
Poisson brackets {. , .}m and {. , .}^. In contrast with the Poisson manifolds we have two 
opportunities for the action of the mapping f : M. ^ M preserving both Poisson brackets 

{■,-}m^ — ^{■i-}m {■^■}m {■^■}j\f 

1. /: , 2. /: 

{■■>-}'m^ — ^{-j-IV {■^■Ym {-^-Ym 

In order to distinguish these cases we will say about the Poisson map and the twisted Poisson 
map at the first and the second cases respectively. 

Let us consider the Euler equations on the Lie algebra so(4) 

s = sx VsH{s, t), i = tx WtH{s, t). (5.23) 
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Here vectors s and t are coordinates on so(4) = so{3) © so{3) with the standard Lie-Poisson 
brackets 

The brackets ()5.24|) respect two Casimir functions 

A* = \s\^, S* = |tp, (5.25) 

which are therefore integrals of motion for (j5.23p in involution with any function on the 
phase space. 

Equations ()5.23|) describe the motion of a rigid body with elliptic cavities filled with 
ideal fluid if the Hamilton function is quadratic form 

H{s, t) = {s, Ais) + {S, A2t) + {t, Agt). 



with the special numerical matrices A^ 

In ^Sj Steklov found the Hamilton function for which equations ()5.23|) possess a fourth 
additional integral 

H = c,Hi + C2H2, |ifi,if2}* = 0, 

where ci_2 are numerical parameters and 

Hi = 2{VA^s,t)-{s,As), H2 = {t,A''t)-2{VA^s,At). (5.26) 

In pP Bobenko established isomorphism between the Steklov system on so(4) and the Steklov- 
Lyapunov system on e(3). 

Proposition 5 /i/ // the phase space M ~ so(4) is identified with the space M ~ e(3) by 
the following linear map 

?■■ (P, M)-.(s = 2p, t= -^{M- Bp)^ , (5.27) 

where A, B are given by / ti.,5pi..9)) . then the equations of motion 115. 2'^) for the Steklov system 
on so(4) coincide with the Kirchhoff equations M.'J^) for the Steklov- Lyapunov system on e{3), 
that is, 

Here we changed p —>■ 2p and M —>■ 2M in comparison with original map pQ to make formulas 
for the Poisson pencil slightly more symmetric. 
We have to underline that 



H^{s,t)=^f(^H,ip,M)y 



as for the usual canonical or Poisson transformations. In fact the map / fj5.27|) gives rise the 
second Poisson brackets on so(4) 

{si,Sj f^ = , [si , tj J2 = Eijk ^—- , [ti , tj I2 = Eijk I ^ ^-—- I , (5.28) 
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whereas inverse map / ^ generates the second Poisson brackets on e(3) 

{Mi, Mj}2 = Sijk ( AfcfcMfc I^Pfc j , {Mi,pj}2 = Sijk—^Pk, {Pi,Pj}2 = ^ijkY ■ 

(5.29) 
Here C^ is the cofactor matrix to matrix C (jl.5j) . 

Proposition 6 The brackets {., .}i il.^) and {., .}2 \5.29^) on the manifold e{3) and brackets 
{., .}l l[5.24}) ^^^ {•! •}2 (E13' on the manifold so(4) are compatible. 

The proof consists of the verification that every hnear combination of these brackets is still 
a Poisson bracket on e(3) and so(4) respectively. 

This Proposition allows us to say that the mapping / ()5.27|) is a twisted Poisson map, 
which identifies two bi-Hamiltonian manifolds e(3) and so(4) such that 

^XT instead of as for the usual Poisson map. 

^2 ^2 ^2 ^^ ^2 

The bi-Hamiltonian systems are defined on the bi-Hamiltonian manifolds by the coefficients 
of the Casimir functions of the Poisson pencil {. , .}a = {• , -ji — A{. , .}2. In our case on the 
bi-Hamiltonian manifold e(3) cubic polynomial PsiX) = |^P ()2.2|) 

3 

P3(A) = J2 ^kX'' = ^-^^ + (^ - ^tr A)\^ + X ^ + T ■ ^^'^^^ 

fc=0 

is a Casimir of the Poisson pencil 

^^ = ^2-X^i, ^xdP3{X,p,M) = 0. (5.31) 

Here ^i and ^2 are the Poisson tensors associated with the brackets {., .}i (jl.2|) and {., .}2 

On the bi-Hamiltonian manifold so(4) polynomial 

3 

P3(A) = /(P3(A)) =J2'^kX'' =A*X^ + HiX^ + H2X-detAB* (5.32) 

fc=0 

is a Casimir of the Poisson pencil 

^* ^ ^* _ x^*^ j^* dP^lX, s,t) = 0. (5.33) 

Here ^i and <^| are the Poisson tensors associated with the brackets {.,.}i (j5.24p and 

As usual, coefficients of the polynomial PsiX) (j5.3Up form a bi-Hamiltonian hierarchy on 
e(3) starting from a Casimir TYq of J3^2 and terminating with a Casimir of <^i 

^2dnQ = 0, ^2dn^+i = J^idHi, ^idHa = . 
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Coefficients of the polynomial i^3(A) (j5.32p form the similar Lenard chain on so(4). Therefore, 
the Steklov-Lyapunov system on e(3) and the Steklov system on so(4) are the Gelfand- 
Zakharevich systems 4j. 

It is known that the Gelfand-Zakharevich systems admit the Stackel separation of vari- 
ables (see jH] and references within). However in our case tensors ^1,2 and ^12 ^^^ degener- 
ate and, therefore, the corresponding separated variables could be obtained after a suitable 
reduction of the Poisson structures E] only. 

A natural way to do that is to fix the values of the Casimir functions of ^\ or !^\. 
However, in the our case each Poisson tensor ^\ (^1 ) or ^2 {^2) can be properly restricted 
to a corresponding symplectic leaf, but the other tensor does not restrict to same leaf. So, 
a quite general reduction technique given by the Marsden-Ratiu theorem El have to be 
applied in this case. 

In fact we make the necessary reduction when we introduce the chart (x, J) related to 
{p, M) by the map fg (I3.1(jp . Namely, composition of the maps / (I5.27p and fg (j3.1Up gives 
rise to the second Poisson brackets on the second sample of e(3) (J3.H) 

\Jii Jj J2 ^ ^ijkO'k'Jki \"-'i) ■^j J2 ^ ^ijk^j^k ^ii j j"^ j \-^ii -^j J2 ^ " (^O.o4J 

The Casimir of the Poisson pencil ^( — X^2 coincides with the cubic polynomial -P3(A) 

PsiX) = SX"" + 7A2 + H,X + H2, ^; dPsiX, x,J)=0. 

This polynomial Casimir includes two numerical constants 7 and 6 instead of the Casimir 
functions and, therefore, the corresponding symplectic foliation is trivial. According to p!H] . 
in this case the separated variables ()4.15j] are the solutions of the following system of algebraic 
equations 

|-P3(A,x,J) = 0, ^Ps{X,x,J) = 0. 

These separated variables have to coincide with the eigenvalues of the reduced Nijhenuis 
tensors ^ = ^\^2^ or jV = S^\S^2^ . As sequence these variables have to be in the 
involution with respect to both Poisson brackets similar to the integrals of motion. We can 
check this property of the Kotter separated variables (|2.7|) without reduction of the Poisson 
tensors. 

Proposition 7 The vectors 

(M -Bp) xp , _ y/A^t X s 

X = TT-z — ;^^ r and X — 



|(M-Bp)xp| VA^txs 

and, therefore, zeroes of the standard equations 

(x, (A - A)"-^a;) = 0, and (x, (A - A)"^x ) = 

are in the involution with respect to both Poisson brackets on e(3) and so(4) respectively. 

This Proposition has a straightforward consequence. The roots M1.2 = f{ui^2) of the non- 
diagonal entry of the corresponding rational Lax matrix ()4.1(ij) 

= (\-^;)(^-f^) = E T^ = (5.35) 

^ ' det A- A ^^ A-tti ^ ' 

^ 1=1 
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are in the involution {ui,U2}i = {^11^2)2 = on so(4) and, according to Proposition El 
they satisfy to the following equations 



d^ 2JF3(Mi)det(A-^i) d^ 2JF3(w2)det(A-«2) , , 

■^ui = ^ > :^^2 = • (5.36) 

at ui — U2 at ui — U2 

Hence Ui^2 are the separated variables for the Steklov system on so(4). These variables are 
poles of the eigenvector ^(A) of the Lax matrix =Sfe(A) 

■^e(A)ij = Eijky X — ttk {sk + A A/A^tfc 
with the following dynamical normalization 



The conjugated momenta t?i_2 are easy defined by the diagonal entry of the corresponding 
rational Lax matrix ()4.16|) 



V, = -^ {H,,e{X)} 



{vj,Uj}l = 6jk 



X=Uj 



after calculation of the relation ()4.19|1 with respect to the second Poisson brackets ()5.34|l . 
The corresponding separated equations are equal to 

y^ + ^^ = 7 = 12 

such that the Stackel matrix 

S=( ""' ""H (5.37) 

\Ui U2 J 

is another block of the same Brill- Noether matrix Uq (I2.13p . According to ^Tj the Stackel 
systems associated with the different blocks of a common Brill-Noether matrix are related 
by canonical transformation of the time. 

In contrast with the Steklov-Lyapunov system for the Stackel matrix ()5.37|) associated 
with the Steklov system we do not know how to introduce coordinates in which equations 
of motion are the Newton equations. 

6 Concluding remarks 

The first result in this paper is that, starting with the Kotter separated variables we construct 
the Poisson map which transforms the Steklov-Lyapunov case of the Kirchhoff equations to 
the Newton equations on the sphere. A natural way to construct such maps is to identify 
the separated variables and the corresponding separated equations. 

The separation of variables for the Steklov-Lyapunov system on e(3) and for the Steklov 
system on so(4) is discussed in framework of the Sklyanin method and in the bi-Hamiltonian 
approach. The main unsolved questions are how to construct a suitable normalization of 
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the Baker- Akhiezer function and a suitable reduction of the degenerate Poisson tensors. To 
solve these questions on this example we could to construct unknown separated variables for 
the Clebsch system on e(3) and the Shottky-Manakov system on so(4). As above, the bi- 
Hamiltonian structure of these manifolds are defined by the similar linear change of variables 
proposed in jTj. 

As a last remark we recall that in [3] there are separated variables, the elliptic and rational 
Lax matrices and the compatible Poisson tensors for the multi-dimensional Steklov systems. 
In our opinion these results deserve further investigation in framework of the separation of 
variables theory and of the general reduction theory for bi-Hamiltonian manifolds. 

The research was partially supported by RFBR grant 02-01-00888. 
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